Phase diagrams of flux lattices with disorder 
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We review the prediction, made in a previous work [Phys. Rev. B 52 (1995)], that the phase 
diagram of type II superconductors consists of a topologically ordered Bragg glass phase at low fields 
undergoing a transition at higher fields into a vortex glass or a liquid. We estimate the position 
of the phase boundary using a Lindemann criterion. We find that the proposed phenomenology is 
compatible with recent experiments on superconductors. 



It is remarkable that after a decade of experimental 
and theoretical efforts, the phase diagram of type II su- 
perconductors in a field is far from being completely elu- 
cidated Stimulated by the discovery of the high Tc 
materials, a reexamination of the mean-field phase dia- 
gram unravelled two main new phenomena. First, it was 
realized [||J^] , and observed B that due to enhanced ther- 
mal fluctuations the Abrikosov lattice melts well below 
H C 2 into a to a flux liquid. Secondly, it was argued that in 
the solid phase, pointlike disorder could produce a glassy 
state with barriers U(j) diverging at small j, and thus 
characterized by the true vanishing of the linear resistiv- 
ity even at finite temperature |&^| . This was a significant 
departure from traditional models of thermally assisted 
flux flow, which assumed finite barriers between pinned 
states. A precursor sign of an instability towards a glass 
was also found in the flux liquid @. Both for the tech- 
nological applications of high-T c materials and from a 
purely theoretical point of vue, the understanding of the 
detailed properties of such a glassy phase is of paramount 
importance. 

Two main phenomenological theories have been put 
forward to describe this glassy phase and to account 
for some of its properties observed in early experiments, 
mainly the observed continuous transition from the glass 
to the liquid and the giant thermal creep. The first ap- 
proach is based on the gauge glass model ||||, and as- 
sumes a complete destruction of the Abrikosov lattice. 
The second approach retains the the elastic lattice struc- 
ture at small scale |J. Although different in nature, both 
theories agreed that the disorder essential to produce the 
glassy low temperature phase and the vanishing of the 
linear resistivity, was also destroying at large scales the 
perfect flux lattice existing in mean field theory. The low 
temperature phase was therefore generally expected to 
be a topologically disordered phase, lacking translational 
order. Several calculations supported this point of view. 



Elastic theory predicted at best a stretched exponential 
decay of translational order j6|,p|JTc|] (i.e. a power law 
grow of displacements) and general arguments tended to 
prove that disorder would always favor the presence of 
dislocations The vortex lattice seemed to be buried 
for good. 
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FIG. 1. The stability region of the Bragg glass phase in the 
magnetic field H, temperature T plane is shown schematically 
(thick solid line) . Upon increasing disorder the region shrinks 
as indicated by the thin solid line (see text) . The melting line 
of the pure system is shown as a dotted line, and the vor- 
tex glass transition line (or crossover to the pinned liquid) is 
shown thick dotted. 



A few points were not naturally fitting in the frame- 
work of these theories. Experimentally, a first order 
transition between the glass phase and the liquid was 
observed at low fields [ pl]p^ ] rather than the predicted 
continuous transition, observed at high fields. Also, dec- 
oration experiments of the flux lattice at very low fields 
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(60 G) in several materials were showing remarkably large 
regions free of dislocations (l3) . On the side of theory, 
old calculations on the related disordered elastic random 
field XY model |Q as well as more recent scaling ar- 
guments for the vortex lattice [[l5| suggested, within a 
purely elastic decription, a slower, logarithmic, growth 
of deformations. 

In a recent work we have obtained the first quantitative 
theory of the elastic vortex lattice |lf|[L7| in presence of 
point disorder. Contrarily to previous approaches, it pro- 
vides a description valid at all scales and demonstrates 
that while disorder produces algebraic growth of displace- 
ments at short length scales, periodicity takes over at 
large scales and results in a decay of translational order 
at most algebraic Jll| . One striking prediction is thus the 
existence of a glass phase with Bragg diffraction peaks ! 
This result was derived within an elastic theory, assuming 
the absence of dislocations. However, the very result of 
our calculation, i.e that quasi long range order survives, 
led us to advocate that dislocations would be much less 
relevant than commonly assumed The alleged im- 

portance of dislocations in a disordered system ]8| made it 
mandatory to further investigate carefully whether dislo- 
cations would indeed modify the above result. The strik- 
ing result JT7[ that we found based on energy arguments 
is that dislocations are not favorable for weak disorder 
in d = 3. This implies self-consistently the existence of 
a thermodynamic glass phase, as far as energy and very 
low current transport properties are concerned, retaining 
a nearly perfect (i.e. algebraic) translational order. Since 
this phase exhibits Bragg peaks very much like a perfect 
lattice it was christened the "Bragg glass". Because it 
retains a "lattice" structure and Bragg peaks, this glass 
phase is radically different from the vortex glass picture 
based on a random gauge model. In particular, since 
such a phase is nearly as good as perfect lattice as far 
as translational order is concerned, it is natural to ex- 
pect it to melt through a first order phase transition. We 
proposed fT?] ] that the phase seen experimentally at low 
fields was in fact the Bragg glass, solving the apparent 
impossibility of a pinned solid. This allowed to account 
naturally for the first-order transition and the decoration 
experiments. Our prediction filf l that a new phase with- 
out topological defects, should be stable at weak disorder 
received subsequent further support both from numerical 
simulations |2ip2| and from analytical calculations in a 
layered geometry |23]p4||. 

Once the existence of a weak disorder/low field Bragg 
glass phase is established the question arises of deter- 
mining its limits of stability and phase boundaries. The 
Bragg glass phase should be stable as a self-consistent so- 
lution in the elastic limit, i.e as long as R a > a p7| . This 
condition is violated when the field is increased and we 
proposed in |ri| that upon raising the field the Bragg 
glass should undergo a transition into another phase, 
which could be a pinned liquid or another glass (vor- 



tex glass). A natural possibility then was that the crit- 
ical point occuring on the melting line |l2| was the end 
point of the transition line between the Bragg glass at 
low fields and a topologically disordered glassy phase (or 
a strongly pinned liquid) at higher field. We pointed out 
that the fact that this point can be lowered by adding 
impurities was a hint that it was related to this transi- 
tion. Such a field-driven transition corresponds to the 
destruction of the Bragg glass by proliferation of topo- 
logical defects upon raising the field, which is equivalent 
to increasing the effective disorder. The other transition 
from the thermal liquid into the putative superconduct- 
ing state at higher fields is presumably continuous. The 
topology of the phase diagram proposed in fl?}] is as de- 
picted in Fig. 1. 

Several recent experiments can be interpreted to con- 
firm the picture proposed in flT7| . Neutron experiments 
can be naturally interpreted in term of the Bragg glass 
p6| , p7t . In BSCCO neutron peaks are observed at low 
fields and disappear upon raising the field [^8| . The phase 
diagram of BSCCO has been recently explored in details 
by overcoming spurious effects due to geometrical bar- 
riers p9| . It can increasingly be interpreted as a con- 
firmation of our theory, the so-called second magnetiza- 
tion peak line being the candidate for the predicted field 
driven transition. Since our proposal this line has been 
investigated in more details in BSCCO |50| and found 
to be relatively temperature independent at lower tem- 
peratures and to be shifted downwards upon increase of 
point disorder through electron irradiation f3l| . Also, 
similar types of phase diagrams are observed in a vari- 
ety of materials, including YBCO The fact that 
a controlled increase of point disorder through electron 
irradiation shifts the transition line to lower fields [^5| is 
a strong indication that our picture is relevant in these 
materials as well. 

In the present paper we follow on the theory exposed in 
]l7t . We make more quantitative estimates of the phase 
diagram depicted in Fig. 1 using a generalized Linde- 
mann criterion. We also explore in more details some 
experimental consequences of our theory. 

Let us consider a vortex lattice system in presence of 
disorder. We can model the vortex lattice by stacks of 
coupled planes. The system is therefore described by lay- 
ers of two dimensional triangular lattices of vortices. We 
denote by Ri the equilibrium position of the vortices in 
the absence of disorder, labeled by an integer i, in the xy 
plane, and by u(Ri, z) their in-plane displacements which 
are two dimensional vectors, (the vortex can only move 
within the plane), z is the coordinate perpendicular to 
the planes and along the magnetic field and x — (r,z). 
The total energy is: 
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c±4:{dzU a ) 2 ] + J d 2 rdzV(r,z)p(r,z) (1) 

where the density of vortex lines is simply defined by 
p(x) = J2 t 5(2) (r-Ri- u(Ri, z)). The last term in (§) is 
the coupling to disorder. In the limit where many weak 
impurities act collectively on a vortex, point disorder can 
be modelled by a gaussian random potential V(x) with 
correlations: V(x)V(x') = A(r — r')8(z — z') where A(r) 
is a short range function pQ | of range £ (the supercon- 
ducting coherence length) 

A(r) = d U 2 e~ r2 /Z 2 (2) 

where d is the distance between layers and U p a typical 
pinning energy per unit length along z. 

In the high-T c Abrikosov lattice, one has in principle 
to use non local elasticity, and the calculation along the 
lines of ji~7| can be done simply by using the complete 
known wavevector dependent expressions of the 

elastic constants in (Q). Since we are only interested in 
nearest neighbor correlations and want to obtain only an 
order of magnitude of the scales involved we use sim- 
ple constant elastic moduli. The physical properties of 
(|l|) were examined in details in Jl6|,[l7j and we just recall 
here the results needed for the phase diagram. The mean 
squared relative displacements of two vortices separated 
by a distance r is: 

B(r) = ((u(0,0) - u(r,0)) 2 ) (3) 

where () denotes the thermal average, whereas — is the 
disorder average. From B(r) one defines two length scales 
R c and R a in the xy plane (and similarly L c and L a along 
z) such that Bdi S (R c ) ~ max(£ 2 , (u 2 )t) (see below) and 
B(R a ) ~ a 2 respectively. R c is the Larkin-Ovchinikov 
pinning length p3| directly related to the critical cur- 
rent, whereas R a is the scale at which one enters the 
asymptotic regime with a logarithmic growth of the dis- 
placements. The model (|l]) leads to the Bragg glass phase 
with quasi long range translational order. 

To determine the region of stability of the Bragg glass 
phase, we follow the arguments proposed in jljj that the 
elastic structure will become unstable when the displace- 
ment between two neighbors becomes of order of the lat- 
tice spacing a, i.e 

B(r = a) ~ a 2 (4) 

To be more quantitative, one can introduce, as for the 
normal thermal melting a Lindcrmann constant cl and 
take for the criterion of stability of the Bragg glass phase 

B(r = a) = ((w(0,0) - u(a,0)) 2 ) = c\a 2 (5) 

cj,, the Lindemann constant is usually of the order of 
cl ~ 0.1 — 0.2 in usual melting and we make here the 



assumption that is indeed a constant independent of 
the field. 

From (^]) one sees that both disorder and thermal fluc- 
tuations act together to increase the displacements. In 
fact the formula (3.18) of shows that B(r = a) 

splits naturally in two parts B(r = a) w 2(u 2 )t + 
Bdis{f = ci)- One immediate consequence of (^|) is there- 
fore that the melting line should be pushed downwards 
in the presence of point like disorder. In fact the Bragg 
glass can disappear in two ways: (i) if the temperature 
is raised, it will melt to a liquid phase (ii) if the field is 
raised, which amounts to vary the effective disorder in 
the system, the system can become so disordered even 
at short length scales that dislocations will appear. (j|) 
gives thus the limit of stability of the BG phase in the 
H — T plane. Although the complete "melting" curve 
can be computed using the formulas for B(r) obtained 
in Jl(|[l7j], such a calculation is tedious and offers little 
insight. We therefore study mainly here the two limits 
of low temperature where the transition is mainly field- 
driven, and for temperature close to the melting curve in 
the absence of disorder. 

If the temperature is close to the pure melting line, 
B(r = a) is dominated by the thermal fluctuations. Since 
for weak disorder R a 3> a the disorder induced displace- 
ments are negligeable at the scale of nearest neighbors 
and one can compute (||) using thermal fluctuations only. 
One then easily recovers the pure melting line 

T m w 4a 3 v / c^cilc| (6) 

Disorder effects will push the melting line slightly down, 
but effects should be negligible at low field for which the 
effective disorder is small enough. Upon increasing the 
field disorder induced displacements will increase, forc- 
ing the transition line defined by (^|) to go down to zero 
temperature at a finite field Hm- The scale at which dis- 
order dominates can easily be obtained by looking at zero 
temperature. To obtain a reliable order of magnitude of 
the "disorder- induced melting" field Hm, it is necessary 
to know the precise B(r) in the presence of disorder and 
not only its asymptotic forms. Fortunately such a calcu- 
lation was performed in |l6| , |i"7| ]. Using formula (4.18) of 
Ref. JitJ one gets 

B(r) = ^b(r/R a ) (7) 

For r — R a one has from JlTj that b w 1 while for r < R a 
one is in the random manifold regime and one can ap- 
proximate B{r) ~ ^ 3 -(r/i? Q ) 1 / 3 . From the solution of |l7|] 
we know that the above formula is quantitatively correct, 
and not only asymptotics. Using (||) one finds that 

a/R a - (ttc l ) 6 (8) 

Using cl = 0.12 gives R a ~ 350a. Thus the transi- 
tion occurs well before the asymptotic regime. We will 
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find that it does occur (e.g in BSCCO) indeed deep into 
the random manifold regime. One also notes that in 
simplified models without intermediate random manifold 
regime (where one directly go form a Larkin regime to 
the asymptotic regime) the above formula would give 
a/R a = (ttcl) 2 ^ 4 ^^. The transition then occurs for 
smaller values of R a /a, in agreement with the results 
found in a special geometry p3p4 |. 

Using (0) and the expression (4.12) of |17| for R a : 



„ 4 3/2 1/2 

za c 66 c 44 



(9) 



as well as cqq — eo/(4a 2 ) and C44 sa ceo/(7 2 a 2 ) (sin- 
gle vortex contribution) with eo = (&o/4:ir\) 2 and c a 
numerical constant [p2|. One gets: 



4tt 3 



(10) 



One thus obtains an expression for the transition field 
Hm naturally expressed in terms of some characteristic 
fields of the system: 



H M (T = 0) 



(ttCl) 4 , e , 4/3 2 /3„l/3 
" V r r I n c2 12 cross 



(11) 



where we have introduced the crossover field H cmss — 
irc<t> /{j 2 d 2 ) with c ~ ln( 7 d/0 §§ and H c2 = $ /£ 2 - 

As numerical estimate of the melting field Hm for 
BiSCCO with H cmss ~ IT, H c2 ~ 100T, [7" p /e = 0.4, 
cj, = .12 gives Hm ~ 400G in good agreement with the 
observed experimental values pl[ . The fact that this 
field is well below the decoupling field validates a poste- 
riori the calculation ( note also that B(r = 0, z = d) is 
still small at the transition). The general shape of the 
phase diagram is in agreement with the one of Figure 
1. Note that some non linear effects, such as screening 
of disorder by thermal fluctuations or by interactions at 
short scales may not be captured directly by the gaussian 
theory of They can be incorporated by a renormal- 
ization of the effective disorder U P (T). Such effects were 
computed in the flux liquid using RG in M and it was 
shown that the pinning length was renormalized upward 
(and thus the effective pinning strength downward) by 
a factor exp((T/T dp ) 3 ) where T dp ~ ([/ 2 <i£ 2 ceo/ 7 2 ) 1/3 is 
the single vortex depinning temperature pj . It would be 
interesting to compute these effects in the solid as well. 
On general grounds that thermal fluctuations can only 
weaken the disorder, one expects an additional curvature 
upward of the Bragg glass instability line Hm(T) when 
T increases beyond 0(Ta p ). 

In Fig.l two main regions can be distinguished: if the 
temperature is high the stability line is nearly indistigu- 
ishable from the melting line of the pure system. This 
regime corresponds to the case where R a (T — 0) 3> a. 
In that case the translational order is only affected at 



distances huge compared to the lattice spacing, and the 
modification compared of a pure lattice is negligeable as 
far as the melting is concerned. This part of the stabil- 
ity line is therefore nearly identical to the melting of a 
pure lattice and one can expect the transition to be first 
order. The Bragg glass melts to a liquid phase, nearly 
insensitive to disorder. 

If the field is increased one will shorten R a (T = 0). 
The disorder itself is now able to make dislocation pro- 
liferate. In particular even at T = disorder destroys 
the Bragg glass. In this range of field and at low T 
the transition line flattens as a function of temperature, 
since it is controlled mainly by the disorder. The phase 
into which the Bragg glass "melts" at low T is relatively 
poorly understood. It is characterized by an absence of 
translational order and of Bragg peaks. There should 
still be some amount of pinning at low T, but whether 
such a phase is a true glass with diverging barriers, sim- 
ilar to the proposed vortex glass of Ref . || , or simply a 
very viscous form of the liquid phase remains controver- 
sial. This phase could also retain hexatic order (hexatic 
glass) since, at least at a naive level, similar arguments 
for survival of hexatic order as for topological order in 
the Bragg glass can be given. If the phase is a true glass 
phase then it should melt thermally to the liquid, on the 
thick-dotted line of Fig. 1. Whether a true vortex glass 
phase exists in untwinned samples is an important still 
open and controversial question 37 which may need to 

Since 



be settled by high sensitivity ]3S ] measurements, 
the low temperature phase is in any case much more con- 
tinuously related to the liquid phase, one can expect now 
the transition to become second order. The Bragg glass 
therefore provides one natural explanation for a change 
of the order of the (thermal) melting transition, as well 
as for the existence of a field-induced transition. 

One should also point out that at very low fields 
(B ~ H c i) where screening is important a > A, a similar 
(inverted) field driven transition should also occur when 
the field is lowered, from the Bragg glass to a pinned 
liquid (or another glass) as suggested by the decoration 
images of jH| . As shown in jjj the liquid becomes unsta- 
ble for B < Boexp(—(T/Tdp) 3 ), presumably to a glass or 
a pinned liquid. 

The proposed field-induced transition between the 
Bragg glass and the putative vortex glass being just char- 
acterized by an injection of dislocations, it is not neces- 
sarily linked to a decoupling beween the layers. As a con- 
sequence one does not expect the critical current along 
z to become zero at the field-induced transition, at least 
for low anisotropy system like YBCO. Of course it is al- 
ways possible that in materials with high anisotropy like 
BSSCO dislocations prefer to appear first between the 
planes and the BG-VG and decoupling transition coin- 
cide. Let us however emphasize that it does not need to 
be so and that we also expect our transition to occur in 
purely isotropic systems. Another argument against the 
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field-induced transtion being a simple thermal decoupling 
transition (3^] , is the fact that such a transition could not 
extend down to zero temperature. In any case measure- 
ments of the critical current perpendicular to the plane, 
in particular in YBCO, should help to separate between 
the two effects. 

The suggestion that there may be two different glass 
phases could seem farfetched. There is a case however, 
mostly of theoretical interest at present, where it should 
happen as a direct consequence of our Bragg glass con- 
siderations. This is for d-dimensional vortex line sys- 
tems with correlated disorder or equivalently in d — 1 
dimensional quantum bosons with disorder. It is rea- 
sonably well established theoretically, numerically and 
experimentally that a Bose glass phase exists for these 
systems in d — 3 (i.e d = 2 + 1 (2 space, 1 time dimen- 
sion) for bosons). It is also believed that this phase lacks 
translational long range order in the plane perpendicular 
to the columns. Indeed since the vortices are localized 
along the columns, one can roughly view the properties 
in the perpendicular plane perpendicular, as those of a 
d — 1 dimensional system with point-like disorder. For 
the d = 3 vortex problem, dislocations are therefore ex- 
pected to appear (as they presumably appear for d = 2 
systems with point like disorder). In higher dimension 
however, this need not be the case. For instance in d = 4 
for vortex systems (d = 3 + 1 for quantum particles) one 
is led, by similar arguments as in p"7| to two distinct local- 
ized phases. For weak disorder no dislocation will appear, 
giving a Bose glass with topological order. This "Bragg- 
Bose glass" phase is the equivalent for columnar defects 
of the Bragg glass one occuring for point-like disorder. 
For stronger disorder, dislocations will destroy the topo- 
logical order perpendicular to the columns, giving back 
the "conventional Bose glass" , i.e the continuation of its 
d = 2+1 version. At the transition between these two dif- 
ferent Bose glass phases, unbinding of dislocations loops 
(cylinders) should occur. An interesting point is that 
in the "conventional Bose glass" these dislocations loops 
will remain pinned to the colummnar defects, and it will 
thus be a true glass. This Bose glass phase would be in 
some sense analogous, in the case of point disorder, to the 
putative vortex glass. However, for pointlike disorder it 
is much less obvious that such a phase exists as a genuine 
thermodynamic phase. Let us note that this problem was 
studied analytically within an elastic theory in ]3^| us- 
ing a variational method and in |36| using RG methods: 
the phase described there is thus the "Bragg-Bose glass" . 
The difference between the two phases should occur only 
at scales larger than the distance between unpaired dis- 
locations. 

Several consequences of our theory could be further 
checked in experiments. First since the Bragg glass 
phase has translational order and the vortex glass has 
not, neutron experiments should observe a destruction 
of the Bragg peaks at the same location [|3 than the 



transition observed by magnetic measurements. Such a 
feature seems to be consistent with the existing exper- 
imental data in BSSCO [^8), but a more detailed ex- 
perimental investigation would be needed to check this 
point in other materials as well. Another clear distinc- 
tion between the two phases should be observed when a 
cycling in current similar to the one of ]26| ] is performed. 
Such cycles taking the system above the critical J c and 
then back to zero, are expected to heal the lattice and 
to expel out of equilibrium dislocations. One can ther- 
fore expect good healing in the Bragg glass phase, as is 
indeed the case, since dislocations can only exist as out 
of equilibrium object. On the other hand the same cycle 
performed in the VG phase should make little difference 
on the neutron diffraction pattern since the equilibrium 
state already contains unpaired dislocations. 
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FIG. 2. I — V characteristics in the Bragg glass phase and 
in the vortex glass (or pinned liquid) are shown schematically 
respectively as solid and dashed lines. One goes from small 
(Ji) to larger (J2) critical currents, but rapidly divergent to 
weakly divergent (or finite) barriers, when increasing the field. 

Finally one expects the barriers to vary very differently 
in the two glass phases. In the Bragg glass phase, elas- 
ticity is strong. Pinning can only be collective and one 
expects therefore weak barriers at short length scales. 
This implies a small critical current. On the other hand, 
since creep can only occur collectively, the barrier should 
grow very rapidly with decreasing current. A simple 
estimate gives U(j) ~ (l/j)^ with /1 ranging between 
fi = 0.7 — 0.8 at intermediate currents to /i = 0.5 at 
small currents [jl?]] while taking the dispersion of elas- 
tic moduli into account leads to higher values of fi in 
the intermediate regime [EJ. On the other hand in the 
VG phase barriers should be significantly larger at short 
length scales since the nearly destroyed lattice has ad- 
ditional effective degrees of freedom, such as free dislo- 
cations, and can thus adapt more easily to the pinning 
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potential. The critical current should therefore increase 
when approaching the field-induced transition. On the 
other hand in the VG phase the barriers should grow 
much more slowly with decreasing current since there is 
no need for collective motion, or even remain finite if the 
phase is simply a crossover from the liquid phase. Some 
estimates of the exponents for the gauge glass model gave 
very small exponents of the order of /i ~ 0.1 — 0.2. One 
can therefore expect I — V characteristics evolving with 
fields like the ones shown in Figure 2. Such a behav- 
ior is in good qualitative agreement of the observations 
of a second peak in relaxation measurements when the 
transition is passed. More refined transport, or relax- 
ation measurement should help deciding on the behavior 
of the barriers. 

It is important to note that the Lindemann criterion 
used here is not a detailed theory of the transition when 
dislocation proliferate, which is not yet available. It 
represents one possible mechanism of instability domi- 
nated by short length scales. It thus provides a reason- 
able upper bound for the instability field Hm since the 
Bragg glass certainly cannot self-consistently survive if 
R a < a. However, it cannot be excluded that because 
of the weakening of translational order at large distances 
in the Bragg glass compared to a real solid, unbound 
dislocations start to appear first at large length scales 
providing a different instability mechanism. In that case 
this additional phase (which may or may not be a true 
glass) would also melt through a first order transition 
with good short distance translational order properties. 

In conclusion, we have examined in some details the 
implications, for the phase diagram of type II supercon- 
ductors, of the existence of a glass phase with transla- 
tional order: the Bragg glass. We have shown that it 
provides a natural interpretation of many of the features 
of the phase diagram of BSCCO and YBCO observed in 
experiments. Namely a crossover from a first order melt- 
ing transition to a continuous transition when the field is 
increased, and the existence of a field induced transition. 
We interpret this last transition as the destruction of the 
Bragg glass phase by spontaneous injection of unbounded 
dislocations, into a topologically disordered glass phase 
or liquid. This transition, being disorder driven, should 
extend down to T = 0. 
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